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PREFACE 

The determination of the albedo of Earth is currently of considerable 
interest to atmospheric scientists, meteorologists, and global climatologists. 

The purpose of this study was to investigate the role of particulates in 
the atmosphere and in what way they contribute to the albedo of the Earth- 
atmosphere system, making use of SKYLAB data acquired over specific test 
sites . 

The SKYLAB data and auxiliary aircraft data were analyzed and compared 
with results of calculations based upon current atmospheric-radiative- 
transfer models. 

This work was performed under Contract NAS9-13279 for the National 
Aeronautics and Space Administration. The author wishes to thank the 
project monitor, Z. Byrns and his assistant, W.J. Johnson for the many 
discussions concerning the handling and analysis of data. 

ERIM personnel who participated in the technical part of the investi- 
gation were Peter Lambeck and Lynn Ziegler. The typing of the report was 
done by Jan Dixon. The work was performed in the Infrared and Optics Di- 
vision under the direction of Richard R. Legault; The Principal Investigator 
was Robert E. Turner. 
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1 

SUMMARY 

In recent years, meteorologists, atmospheric scientists, and applied 
mathematicians, aided by modern high-speed computers have been able to 
improve upon the development of atmospheric models. These models of the 
atmosphere are used to predict future weather patterns over an extended 
area and they are also used to study the climate over major portions of 
the Earth for many centuries in the past and on into the future. This 
latter study, global climatology is more than of academic interest, for it 
has much to say concerning the quality of life for all mankind. 

Besides the ravages of nature such as tornadoes, hurricanes, drought, 
floods, and other weather-related phenomena, the seemingly subtle changes in 
the composition of the Earth’s atmosphere can have a profound impact upon 
our existence. The impetus for the present investigation is the concern 
over the long-term change in the composition of the semi -permanent atmospheric 
component called the aerosol. 

The first part of this report deals with general aspects of radiation 
balance for the Earth and the development of atmospheric-radiative-transfer 
models used in the description of the Earth's natural radiation field. The 
last part of the report deals with the calculations of the Earth's albedo 
using the models of the atmosphere and the multispectral data obtained by 
SKYLAB. 
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2 

INTRODUCTION 


The temperature of a planet is determined by equating the incoming flux 
from the sun (multiplied by the fraction of energy absorbed) to the outward 
radiated flux. The temperature is thus dependent upon the solar flux and 
the albedo of the planet, where we can define the albedo as the ratio of 
the reflected radiation to the incident radiation. lit this sense, the albedo 
of the combined atmosphere-surface system is somewhere between zero and 
unity. If the albedo is zero, the planet absorbs all incident radiation 
and re-radiates it as thermal radiation. On the other hand, if the albedo 
of a planet is unity, then none of the incident radiation is absorbed and no 
temperature can be stated. 

The radiation balance of a planet is determined by the solar flux value, 
the reflectance of the planetary surface, and by the reflectance and absorption 
properties of the planetary atmosphere. As a result of man's activities such 
as the large-scale cultivation of the soil, irrigation, and the building of 
roads and cities, the surface features of Earth have been altered over the 
years. This has the effect of changing the surface reflectance or albedo. 
Likewise, man has succeeded in changing the atmosphere by introducing artificial 
gases and particulates which can change the atmospheric albedo. Although man 
has no direct influence on the solar flux, it is known that the energy radi- 
ated by the sun is not constant in time. All of these effects can change our 
environment, either in a direct way by allowing a greater or lesser amount 
of solar radiation to reach the Earth's surface, or in an indirect way by 
altering those atmospheric conditions which lead to changes in our climate. 

It is therefore important that we know how rapidly and in what way the 
planetary albedo of the Earth is changing so that climatologists can improve 
the predictive capability of their models. 

In this report we shall primarily be concerned with the possible vari- 
ation in the particulate component of the atmosphere. Using radiative- transfer 
theory, an atmospheric radiation model has been developed which is then used 
to calculate the spectral albedo of the Earth for a great variety of possible 
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atmospheric states. Interesting relationships are shown which connect 
the albedo to easily measurable quantities for various atmospheres and 
energy conservation relations are defined for any atmospheric state. Finally, 
using SKYLAB multispectral data for specific target areas we attempt to 
describe the aerosol contribution to the Earth's albedo. In doing this we 
clearly point out the difficulties encountered in this method of albedo 
determination and describe techniques for overcoming such problems in the 
future. 
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3 

THE ENERGY CYCLE 

With the exception of radioactive materials inside the Earth, almost 
all of the external energy we have comes ultimately from the sun. Energy 
from the sun is reflected and absorbed by the atmosphere and is reflected 
and absorbed by the Earth’s surface. In this section we shall consider 
the spectral character of the solar radiation and the terrestrial radiation 
within realistic atmospheres. 

3.1 Solar Radiation 

There have been many measurements of the so-called solar constant, i.e., 

the total energy incident upon the atmosphere integrated over the entire 

spectrum. Some of the many references to this work over the years is given 

[ 1 ] 

by Drummond and an excellent history of the measurement of solar radiation 

[ 2 ] 

is given by Henderson . Most of these investigations however, involving an 
extrapolation to zero air mass are considered to be unreliable. More recent 
measurements have been made aboard balloons, high-altitude aircraft and 
spacecraft by Arvesen et al. , and Thekaekara and Drummond^^ which 
indicate that the earlier values of the solar spectrum were too high. 

There could be intrinsic fluctuations in the solar spectrum other than 
the well known large fluctuations in the ultraviolet part of the spectrum. 

We shall use the NASA standard values given by Thekaekara and Drummond for 
our analysis. A portion of the spectrum in the near-UV, visible, and 
near-IR is illustrated in Figure 1. Figure 2 depicts the spectrum 
according to Pettit ^ over an extended part of the spectrum and it also 
shows the spectrum at sea level. Nimbus F satellite, launched in the 
summer of 1975 is providing data on the extra-terrestrial solar spectrum 
in connection with the Earth Radiation Budget investigation and Nimbus G 
will also monitor the solar spectrum. As a result of these future satellite 
measurements , we hopefully will have a more detailed knowledge of the solar 
spectrum and its variations. 
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3.2 Atmospheric Structure 

In this section, we shall be concerned with the basic structure of 
the atmosphere insofar as it applies to the problem of the transfer of 
radiation. 

3.2.1 Atmospheric Gases 

The Earth’s atmosphere consists primarily of molecular oxygen and 
nitrogen and the atomic gas argon. Variable gaseous components are water 
vapor, carbon dioxide, ozone, sulfur dioxide, hydrogen sulfite, hydrogen, 
and other complex molecules such as nitrogen compounds, methane, and 
formaldehyde. The strongest absorption of radiation by gases takes place 
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in the ultraviolet and infrared part of the spectrum. Oxygen (C^) , 
nitrogen (H^) , and ozone (0^) absorb primarily in the ultraviolet whereas 
water vapor (E^O) > carbon dioxide (CC^) and ozone absorb strongly in the 
infrared. Throughout the visible part of the spectrum, the only signifi- 
cant absorption by gases occurs in Chappuis band of ozone near 0.6pm. 

Besides absorbing radiation, gases and small fluctuations in the 
density of air can also scatter radiation. If the scattering center is 
much smaller than the wavelength of the incident radiation, then an oscil- 
lating electric dipole field is created. Lord Rayleigh calculated this 
effect almost a century ago and found the angular and spectral dependencies 
of the radiation to be as follows: 


Intensity « 


X ^(l+cos^G) 


( 1 ) 


where X is the wavelength and 0 is the angle between the incoming and out- 
going radiation. This type of scattering is referred to as Rayleigh 
scattering and is responsible for the blue color of a clear sky. There are, 
in addition, other scattering processes associated with the interaction of 
higher energy photons with gases but these are not of interest to us here. 

3.2.2 Aerosols 

Small particles arising from a great variety of natural and anthro- 
pogenic sources can temporarily be suspensed in the atmosphere. This 
se mi —permanent suspension of liquid or solid particles in the atmosphere 
is called an aerosol. Examples are hazes, clouds, fogs, mists, smoke, smog, 
and dust. The particles form around a tiny nucleus composed of various 
organic and inorganic materials such as ammonia, nitrates, sulfates, and 
hydrocarbons. The sources of aerosols have been summarized by Hidy and 
Brock ^ and are presented in Table I. 

[7] 

The sizes of particles can vary over a wide range. Junge has 
divided the aerosols into three categories j Aitken nuclei with radii 


REPRODUCBXLOT OF TUB 
ORIGINAL EMH1 IS POOR 



FORMERLY WILLOW RUN LABORATORIES. THE UNIVERSITY OF MICHIGAN 



TABLE I. SOME IMPORTANT SOURCES OF ATMOSPHERIC AEROSOLS. 

Production Rate in Tons Day~l on a Worldwide Basis 


(a) Natural 



Estimated 

Max. % 


Estimated 


production 

by wt. of 


production 

Source 

rate 

total 

Source 

rate 

Primary 

2 X 10 4 - 10 5 

2 . 

Secondary 


Dust Rise by Wind 

9.3 

Vegetation 

5 x 10°- 

Sea Spray 

3 x 10 6 

28. 

(hydrocarbons- 

terpenes) 

3 X 106 

Extraterrestrial 
(meteoritic dust) 

50-550 


Sulfur Cycle 
(oxidation of 

10 5 - 10 6 

Volcanic Dust 

10 4 

0.09 

H 2 S— SO 42 -) 


(intermittent) 



Nitrogen Cycle 


Forest Fires 

4 x 10° 

3.8 

Ammonia 

7 X 10 5 

(intermittent) 



NO x — NO 3 

10 6 




Volcanoes (volatiles, 

~ 10 3 




SO 2 and H 2 S) — inter- 
mittent 





Sub-total 

10.1 x 10 6 


(b) Anthropogenic 


Source 

Estimated 

Production 

rate 

Max. % 
by wt. of 
total 

1. Primary 

Combustion and industrial 
Dust rise by cultivation (inter- 
mittent) (U.S. only) 

1-3 x 10 5 
10 2 -103 

2.8 

0.009 

2. Secondary 

Hydrocarbon vapors (incomplete 
combustion, etc. 

7 x 10 3 

0.065 

Sulfates (oxidation of S0 2 and 

h 2 s) 

3 x 10 5 

2.8 

Nitrates (oxidation of NO x ) 

6 X 10 4 

0.56 

Ammonia 

3 x 10 3 

0.028 

Sub-total anthropogenic 

6.7 X 10 5 

7.1 - 6 % 

Total: All sources 

10.7 x 10 3 tons day " 1 


Max. °/ ( 
by wt. (■ 
total 


28. 

9.3 


6.5 

8.3 

0.009 

94. 
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-7 -5 -5 

between 10 cm and 10 cm, large particles with radii between 10 cm 

-4 -4 

and 10 cm, and giant particles with radii greater than 10 cm 
roi 

Deirmendjian L has considered various mathematical formulas for re- 
presenting the size distributions. One of these formulas is the so-called 
modified gamma distribution given by 


n(r) 


6 Y 

ar exp(-br ) ; 0 £ r < 


( 2 ) 


where the various parameters of the distribution function are related to 
the number density and the mode radius of the particles. In Equation (2) 
n(r) is the number of particles per unit volumes of radius r. Figure 3 
illustrates the size distribution functions for haze M (marine or coastal) , 
haze L (continental) , and haze H (stratospheric) . 

For atmospheric optics, the composition of particulates manifests 
itself in terms of the complex index of refraction. It is represented 
as follows: 


m(A) = m^(A)-im^(X) 


(3) 


where m^(A) is the real part of the refractive index which is responsible 
for the scattering of the radiation and m^(A) is the imaginary part which 
corresponds to absorption of radiation. To determine the refractive index 
for an extended region of the atmosphere is not an easy task. Usually in situ 
sampling is performed and absorption methods are carried out in the labora- 
tory to find the imaginary part. In general, as the amount of water vapor 
increases the particles have more water in them and the index therefore ap- 
proaches that corresponding to pure water. Thus, in the visible region 

m(A)-KL. 33 as the water vapor increases. A more detailed treatment of 

rgi 

aerosol absorption is given by Turner for the case of realistic atmos- 
pheres. 

It is usually assumed that the aerosol particles are spherical in 
shape. This assumption is probably valid for water or liquid particles but 
dry particles can have any shape. It is not clear how the non- sphericity 
affects the optical properties of the particles for a polydispersion, i.e. 
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a large collection of particles with a spectrum of sizes. It is commonly 
assumed that if the particles are of random orientation then the gross 
optical property is. similar to that for a collection of spherical particles. 

3.2.3 Attenuation of Radiation 

The attenuation of an infinitely narrow monochromatic beam of 
radiation can be expressed by the following equation: 

x 

I-^(x) = 1^(0) exp [- / K^(x’ )dx* ] (4) 

o 

where I-^(x) is the intensity of the beam at distance x, 1-^(0) is the 
source intensity, and K-^(x) is called the volume extinction coefficient. 

The extinction coefficient K^(x) is a measure of the total loss of energy 
by both absorption and scattering processes along the path. Thus, we 
can write 


K^(x) = ot^(x) + B^(x) (5) 

where a-^(x) is the volume absorption coefficient and $^(x) is the volume 
scattering coefficient. The extinction coefficient can also be divided 
into components representative of the attenuation properties of the 
medium. Thus , 


K X (x) = k ra < x) + k aa ( x) (6) 

where (x) is the extinction coefficient due to the gaseous component 
RA 

and K.,(x) is the extinction coefficient due to the aerosol component. 

iiA 

Since we will assume in this report that we are dealing with a horizontally 
homogeneous atmosphere we shall consider the variability of the atmosphere 
parameters in the z (vertical) direction only. 
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The coefficients for absorption, scattering, and extinction can be 
expressed as 


a-^(z) = N(z)ci a (A) 

(7) 

3^(z) = N(z)cr s (A) 

(8) 

k x (z) = N(z)CJ fc (A) 

(9) 


where N(z) is the number of independent centers of interaction and cr (A) > 

(A) , and cr^A) are the absorption, scattering, and total cross sections 
respectively. For a pure, aerosol-free atmosphere there is only the gaseous 
part of the attenuation to consider and the absorption cross section is cal- 
culated or measured for whatever spectral region is considered. The Rayleigh 
(gaseous) scattering cross section is given by 


a s (X) = 


8iT 3 (m(A) 2 -l) 2 
3 A 4 N 2 


6 + 3p 
.6 - 7p. 


( 10 ) 


where p is the depolarization parameter (p = 0.035), N is the number of 
molecules per unit volume, and m(x) is the refractive index of air. The 
index m(X) is given by 


m(A) = 1.000064328 + 


0.02949810 
146 - 1 1\ 2 


0.00025540 
+ 2 
41 - 1/A 


( 11 ) 


for dry air at standard temperature and pressure of 15°C and 1013.250 mb. 
The wavelength X in Eq. (11) has the units of micrometers. For any tem- 
perature t and pressure p other than the standard values t and p g , the 
refractive index is given by 

t 

m(X) - 1 + (m (X)-l) ( 12 ) 
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Hence, we can determine the Rayleigh cross section and the Rayleigh scat- 
tering coefficient for any altitude. 

The corresponding cross sections for particulates involve much greater 
complications for several reasons: (1) uncertainty concerning the structure, 

(2) composition, (3) sizes, (4) and shape of the particles. If we make 
the assumption, however, that the particles can be represented by uniform, 
homogeneous spheres, then we can calculate the cross sections by using 
classical electromagnetic theory. These calculations were first performed 
by Lorenz 2.890 and by Mie^^ in 1908 and the resulting formulation 

is known as Mie theory. The scattering, absorption, and total cross section 
for particles of radius r are given by the following: 


°sX (r) 


7tr 2 Q 


2TTr 

2 

x 


2 00 


A=1 


( 2 & + 1 ) 



(13) 


a aX (r) = 7Tr 2 Q a = a^r) - a^r) (14) 

2 00 

a fcX (r) = 7rr 2 Q t = l (2A + l)R e (a + b ) (15) 

x &=1 

where the parameters a^ and b^ are given in terms of Ricatti-Bessel functions, 
and the dimensionless parameter x = 2TTr/X. The factors Q , Q , and Q are 

S3. L 

referred to as the efficiency factors. As the particle size increases 

Q and Q approach the value of 2 whereas Q approaches zero. Hence, very 
s t a 

large particles have cross sections which are twice their geometrical values. 

It is beyond the scope of this report to go into greater detail concerning 

the extensive subject of single-particle scattering theory. A full treatment 

ri21 [131 [141 

of this work is given by Stratton , van de Hulst , and Kerker 

From Eq. (10) it is seen that the spectral dependence of a Rayleigh 

atmosphere is approximately \ ^ . For a hazy atmosphere the spectral dependence 

. -1.3 

is a more slowly varying function of wavelength and is given by ’ ,. 
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For aerosols with any size distribution n(r) the formulas for the 
absorption, scattering, and extinction coefficients are given by 


a, 


^ r\ 


(r)a (r)dr 


=f n(r)a s (r)dr 


*' o 


(r)a t ( r )dr 


(16) 

(17) 

(18) 


Another concept which is useful in radiative-transfer analysis, is 
the optical depth and optical thickness of the atmosphere. The optical 
depth T ^(h) at some altitude h in the atmosphere can be defined as 

>•00 

T^(h) =J K-y (z) dz . 


- A VU, - 

Thus, optical depth, a dimensionless quantity, represents the amount of 
material above altitude h. If we are at the bottom of the atmosphere, 
the optical thickness is then given by 


(19) 


oA 


f 


K^(z)dz . 


( 20 ) 


[15] 


For practical considerations, Koschmieder developed a formula 

relating visual range at sea level to the volume extinction coefficient K 
at the peak of the human eye response at a wavelength of 0.55pm. This 
relationship is depicted in Fig. 4. 

Taking values of the molecular number density for the U.S. Standard 
Atmosphere and the Rayleigh cross section of Eq. (10) we have integrated 
the volume extinction coefficient over altitude to determine the optical 

n 6 1 

thickness of the atmosphere . This is illustrated in Fig. 5. Making 

use of Koschmieder *s formula relating extinction coefficient and horizontal 
visual range and using results of measurements of light scattering in the 
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atmosphere, Elterman ^ ^ has formulated a model for hazy atmospheres. 
This model relates optical thickness, altitude, wavelength, and visual 
range. We have plotted his results in Fig. 6. 


3. 3 Climatology 

If we assume radiative equilibrium for the Sun-Earth-Atmosphere sys- 
tem the equivalent ’’brightness" temperature T can be defined as follows: 

s 


T 

s 


(I-A)Eq 

4cr 


1/4 


( 21 ) 


where E is the extra-terrestrial irradiance (flux) incident at the top of 

the atmosphere, 0 is the Stef an-Boltzmann constant, and A is the albedo of 

-2 

the Earth-atmosphere system. Using a value of 1353 Wm for the solar 
flux we can calculate the temperature. This is illustrated in Fig .7, It 
should be noted that the largest change in temperature occurs when the 
albedo is quite high (A > 0.9). Estimates of the planetary albedo usually 
range from 0.30 to 0.40. Thus, a more detailed plot of the temperature for 
a realistic albedo range is presented in Fig. 8. If we have a planetary 
albedo of about 0.33 then according to Eq. (21) or Fig. 8 we have a tempera- 
ture of 251. 5 °K, a value close to the satellite observations of Raschke 
and Bandeen^ 18 ^. This temperature is, however, much lower than the mean 
surface temperature of Earth which is about 286 °K. The reason is due to the 
effect of the absorption of outgoing infrared radiation, primarily by the 
gases, carbon dioxide and water vapor and by clouds which re-emit radiation 
downward to keep the surface temperature high. This is usually referred to 
as the "Greenhouse" effect. Thus, the Earth's temperature is actually a 
function of the albedo of the Earth-atmosphere system, the incoming solar 
flux, and the composition of the atmosphere. The Earth's surface and hence 
the surface albedo has undergone changes over the years. The building of 
roads and highways alone has taken about one percent of the land area in 
the United States. Large-scale land cultivation and irrigation has also 
changed the surface albedo. In addition, a decrease in the forest areas 
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FIGURE 6. AEROSOL OPTICAL DEPTH t a VERSUS VISUAL RANGE V. Parameter is 
spectral wavelength, A • 
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FIGURE 7 

HYPOTHETICAL EARTH 
TEMPERATURE AS A 
FUNCTION OF ALBEDO 
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could produce serious climatic effects since the amount of carbon dioxide 

in the atmosphere is dependent upon the biomass. 

The role of particulates in determining the temperature of the Earth 

is now being investigated. In fact, the increase in the atmospheric 

particulate load could reflect solar radiation to such an extent that the 

mean temperature of Earth will decrease by 3.5°K over a period of several 

[191 

years. Rasool and Schneider have performed calculations using a two- 

stream model for the transfer of radiation through the atmosphere. The 

predicted decrease in temperature as a function of atmospheric optical 

thickness is illustrated in Fig. 9. Other recent calculations^ 2 ^ also 

indicate that a net heating effect could result from an increase of the 

[ 21 ] 

particulates and also from an increase in the chlorof luorocarbons 
or Freons. Clearly, the problem is far from being solved; more reliable 
data are needed on the question of the increase of those gaseous components 
which enhance the Greenhouse effect and, also better atmospheric models 
should be used in the analysis of temperature changes. 
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4 

THE ATMOSPHERIC MODEL 


In this chapter we shall consider the basic atmospheric-radiative— 
transfer model which was used for the analysis of SKYLAB data and for 
aerosol albedo studies. 

4.1 Radiometry 

In this section we shall define the basic radiometric terms which 
are used in radiative- transfer theory. 


4.1.1 Radiance 

The radiance L is defined as the radiant power that leaves a sur- 
face per unit solid angle and unit projected area of that surface. This 
can be expressed mathematically as 


_ 3 2 P 

3A 3fi cos 0 


( 22 ) 


where P is the radiant power, A is the area, Q, is the solid angle, and q 

is the angle between the direction of the outgoing radiation and the surface 

normal vector. Usually we shall be concerned with the spectral radiance; 

2 

the units will be mW/cm -sr-pm. It is well known that the apparent brightness 

of an object is invariant with distance. The total flux decreases as one 

gets farther from an illuminated source but the amount of power density 
1 

per unit solid angle is constant. Thus the radiance-invariance law which 
for two media with different indices of refraction can be expressed as 


(23) 


n_ 


n. 
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4.1.2 Irradiance 

The Irradiance E is defined as the radiant power incident upon a 
surface per unit area, i.e. 



E - ^ 
E ~ 3A 


( 24 ) 


a quantity which does depend upon distance. We shall use the units of 
2 

mW/cm -ym. for spectral irradiance. For a more complete definition of 

these radiometric terms and their conceptual foundations the work of 
[23] 

Preisendorfer is recommended. 

From Eqs (22) and (23) we see that the irradiance can be found if 
the radiance is known over a hemisphere, i.e. 


E = / L cos 0 dfi (25) 

n 

where is one hemisphere or 2tt steradians. 

4.1.3 Transmittance 

The transmittance of a beam of radiation propagating through the 
atmosphere lies somewhere between zero and unity. For zero transmittance 
no direct radiation passes through the medium and for a value of one all 
the radiation passes through undiminished in intensity. We shall neglect 
gaseous absorption in this report because we are primarily interested in 
the scattering of radiation by gases and aerosols. 

If an infinitesimal, narrow beam of monochromatic radiation propagates 
through a medium the intensity at distance x from the source is 

I^(x) = I^(o)T^(x) (26) 

where I. (o) is the source intensity at wavelength A and T^ (x) is the 
A A 

spectral transmittance at distance x. If the loss of radiation from the 
beam is total, i.e., complete absorption and/or scattering out of the beam 
then the fractional change in intensity is given by the following: 
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A I, (x) 

T^0 = - K X (X) Ax (27) 


where is the probability per unit length for a net loss of radiation 

from the beam and Ax is the distance over which the fractional change 
occurs. As the length Ax approaches zero we have 


dl, (x) 

- s — ■ - K A (x) hW 


(28) 


-f X K x ( x ')dx’ 

I x (x) = I (o) e Jo 


(29) 


Thus , we see that 


the spectral transmittance T. (x) is given by 

A 





K x (x' )dx' 


(30) 


Hence, one of the problems is to determine the probability per unit length 
K^(x) . It can be determined experimentally or theoretically by applying 
the Mie scattering theory for particles. The quantity in the exponential 
in Eq. (30) is called the optical depth of the medium. For a vertical 
distribution through Earth's atmosphere the spectral transmittance would 
be given by 


T 


X 


/* m 

“Jo K x^ z '^ dz ’ 


(31) 


The transmittance versus optical depth is depicted in Fig. 10. 
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4.2 Reflectance 
4.2.1 Non-Lambertian 

In general, radiation incident on a surface may be reflected at 
any angle. If the incident radiation is represented by the zenith angle 
0’( = cos \l') and an azimuthal angle cj) 1 and the outgoing reflected radiation 
is represented by corresponding angles 0 and (j), then the bi-directional 
reflectance can be defined by the relation 


L(y,c|>) 


r 2tt ri 


o Jo 


,4> T ) L (-y' , (j) ' ) dy 1 dcj) 1 


(32) 


Hence, the reflectance properties of a surface are in general represented 
by four independent angles. For specular surfaces there is a known relation- 
ship connecting the incoming and outgoing radiation, i.e. 


y = y 5 ; 4> - 4» f + tt (33) 

or 

P spec (V*,<| ) ,W , ,<|)’) S(y'-y)S(<j>*-cMT) (34) 

Thus, Eq. (32) becomes 

L(y,())) = L(-y, <j> + IT) (35) 

Many surfaces have been studied to determine their bidirectional reflec- 

r 23 1 

tance L , a large collection of which exhibit some specular characteristics. 
4.2.2 Lambertian 

If the surface has no dependence on angles then we say that the 
surface is Lambertian or perfectly diffuse. In this case, Eq. (32) be- 
comes 

L(y,(j>) = p’ E (36) 
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where E is the irradiance on the surface. If we now integrate Eq . (36) 
over the outgoing angles, 0 and <f> , we have 


2tt /»1 


y L(y,4>)dvid4> = 


2TT 


O J 


p' y e dyd4> 


(37) 


or 

M = TTp ' E (38) 

where M is defined as the radiant exitance. By definition the ratio of 
exitance to irradiance is the albedo of the surface. Designating the 
albedo or perfectly diffuse reflectance by p we have 

M = pE (39) 


so that Eq. (32) becomes 


L(y,<f>) = ~ E (40) 

for a Lambertian surface. Although such a surface is unrealistic, the 

approximation is used to simplify analysis. Also, for many natural materials 

such as soil and plants the approximation is probably not too bad 

for the visible part of the spectrum. A compilation of reflectance has been 

r 241 

produced by Leeman et al. for over 3800 materials. 

4 . 3 Radiative-Transfer Theory 

It is beyond the scope of this report to go into the details of 
radiative-transfer theory. Instead, we shall consider various solutions 
of the radiative-transfer equation in order to apply the theory to the 
investigation of the albedo of the Earth-atmosphere system. A more complete 
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[251 

treatment of radiative-transfer theory is given by Anding et al. 
Malila et al. , Turner and Spencer ^ 2 ^, Horvath et al. , 
Turner t 28 # 29 * 30 ^ Turner et al. *- 31 \ and LaRocca and Turner*- 32 ^ 1 . 


4.3.1 The Radiative-Transfer Equation 


Since we are concerned with the steady-state condition of 
the transfer of energy in the atmosphere, we can use the time-independent 
radiative- transfer equation 


A 

n 


• V L(r,fi) + K(r)L(r,n) - J 

6 ’ 


p(r )L(r,fi' )dfi’ 


Q(r,J2) 


(41) 


/\ ^ 
in which is the incoming radiance vector, and is the outgoing radiance 

vector. Since we are not considering polarization, the spectral radiance 
L(r,f2) is a scalar. The quantity K(r) is the volume extinction coefficient, 
j3(r) is the volume scattering coefficient, and Q(r,fi) is the source term. 

The integral term represents the contribution of radiation as a result of 
multiple scattering at point r in the medium according to a particular^angular 
distribution disignated by the single-scattering phase function p(r,£2*£J’). 

The three-dimensional radiative-transfer equation does have an exact solution 
for certain cases L J , , but for our purposes we shall consider a one- 

dimensional equation in which the only spatial variability occurs in the ver- 
tical (z) direction. For solar radiation entering the top of Earth's 
atmosphere with a zenith angle 0 (= cos H* q ) and an azimuthal $ o , the one- 
dimensional radiative-transfer equation is 


to (t) r2n cl 

H = L(t,T1,<|>) - -°- 4 — - j I p(T,]i,()),V f ,cj>’) L (t,U’ ,<j>')du’d(j>' 
to (t) 


(42) 
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in which U) (r) is called the single-scattering albedo and E (t) is the 
o s 

solar irradiance at optical depth T. We are assuming no thermal emission. 


4.3.2 Formal Solutions 

If we assume a vertically homogeneous plane-parallel atmosphere, 
i.e. one in which neither the single-scattering albedo, W (!) nor the 
single-scattering phase function, p (x ;P ,V* ’ f ) depend upon optical depth, 
then the formal solutions of Eq. (42) are the following: 


- (T -T)/y 

l(t, y,4>) = L(T o ,y,<|>)e 


+ 


to 


4TTy 


2tt 1 


’ 

f 

p(y » 4 >,y* , 4 > f ) 

0 

i — 1 
! 


r( T ' T) /y L(T T ,y* ,4) « ) dT: ’ dP ' d4> * 


+ 


Wo p(M *"VV 


4ir(u+y o ) 


-x/y -t /y ~(x -T)/y 
o o o o 
e — e e 


(43) 


for the upwelling radiation, and 


L(x,-y,cf>) = L(0,-y,(j))e 

2 1 


-x/y 


to 


+ 


4Try 


0 


p(-y,4>,y f 


-i 


WoPt-V’t'-Vo’V 

4TT(y-y Q ) 


for the downwelling radiation. 


x 

e -( T _T ' ) /y L (x ' ,y ’ ,(f) f ) dx ' d y McJ) ’ 

0 


-x/y -x/y Q 

e - e 


V 7* y (44) 

o 
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The boundary conditions are 

L(0,-y,<j>) = 0 


(45) 


i.e. it is assumed that there is no diffuse radiation entering the top 
of the atmosphere, and 



r 2TT 

rl f 

l(t q , y,c|>) = 


TJ*p(y,«l>»u f ,4>») i 

« 

0 

0 




+ L SKY (T o’" y ’ ’^J (46) 

Hence, Eqs. (43) and (44) together with the boundary conditions, Eqs. (45) 
and (46) should allow us to determine the complete radiation field in 
the atmosphere. It is interesting to note that Eq. (43) can be written 
in the compact form 


L = L T + L 
o p 


(47) 


where L is the surface (target) radiance, T is the transmittance between 
o 

the surface (target) and the sensor, and L is the so-called path radiance, 
i.e. that radiance which arises from singly or multiply scattered radiation 
into the instantaneous field of view. Equation (47) can also be thought of 
as the remote sensing equation. Likewise, if we are at the bottom of the 
atmosphere (t=t )' then the sky radiance is given by Eq. (44) with t=T q . 


4.3.3 Approximate Solutions 

The radiative-transfer equation, either in its differential form, 
Eq . (42) or its integral forms, Eqs. (43) ,(44) , is quite difficult to 
solve. Exact solutions do exist for special cases ^ 34 ’ 35 , but for the 
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general cases of a realistic atmosphere with an aerosol distribution 
elaborate mathematical techniques must usually be employed to find a 
solution. One of the main reasons why the equation is so difficult to 
solve is that the scattering of radiation by an aerosol particle is highly 
anisotropic; the scattering is predominantly in the forward direction. 

This is illustrated in Fig. 11. Mathematically, the highly anisotropic 
phase function can be represented by a series of Legendre polynominals 
but the number of terms required to represent the function adequately is 
usually quite large. This means that various analytic solutions are com- 
plicated or require a considerable expenditure of computer time. 


Isotropic and Rayleigh Functions 
■ are Multiplied by Ten 
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FIGURE 11. ANGULAR DEPENDENCE OF SINGLE-SCATTERING PHASE FUNCTIONS IN 
ANY AZIMUTHAL PLANE. 
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For these reasons, we shall consider a simpler solution of the 
radiative-transfer equation. For a first approximation let us represent 
the single-scattering phase function by two-streams, one in the forward 
direction and one in the backward direction. Thus, 

p (y , 4) ,y ',<()') = 4*nrsS (y— y ' )6 (4>— <f> * ) 

+ 4tt(1-t">) 6 (y+M ' )<$ (7T-Hj)— cj> T ) (48) 

where r) represents the fraction of radiation which is scattered into the 
forward hemisphere. For Rayleigh scattering n = 0.5, whereas for a hazy 
atmosphere T) == 0.96. If we also represent the radiance by the same set 
of functions, i.e. 


l(t ,y,4>) = 


— E ' (t ) 6 <y — u ) 6 (tt-h(> o — ) + E_'(t)6(y+y o )6(4>-4> 0 ) 
o L 


where E ' (t) and E ' (t) are the upward and downward diffuse fluxes for an 
atmosphere bounded by a black surface. Proceeding in a similar manner 
for non-black surfaces we can then insert expressions (48) and (49) into 
Eq. (42) in order to solve for the fluxes in the realistic atmosphere. 

The solutions are too lengthy to present here; they are given in Appendix I. 
These solutions will be used in the analysis of the SKYLAB data. 

Using the flux equations along with an expression similar to that of 
Eq. (49) then allows us to calculate the path and sky radiances. These 
results are more complicated than those for the fluxes and they are given 
in Appendix II. 
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The accuracy of this model can be tested by comparing the results 
of the model calculations with other mathematically exact calculations and 
with experimental data. An example of a comparison with Chandrasekhar's 
exact calculations is illustrated in Fig 12. A comparison with experi- 
mental values of sky radiance are depicted in Fig. 13. In either case, 
the agreement is very good considering the fact that the model was not 
expected to conform to atmospheres with such a small aerosol content. 



SCAN ANGLE (°) 


FIGURE 12. DEPENDENCE OF SKY RADIANCE ON SCAN ANGLE (PERPENDICULAR TO 

SOLAR PLANE) FOR A RAYLEIGH ATMOSPHERE. Wavelength = 0.546^m; 
solar zenith angle = 36.9°. 


41 




SKY RADIANCE (mW/cm -sr-/im) 


FORMERLY WILLOW RUN LABORATORIES. THE UNIVERSITY OF MICHIGAN 



•80 -60 -40 -20 0 


40 60f 80 

Sun 

ZENITH SCAN ANGLE (deg) 


FIGURE 13. 


DEPENDENCE OF SKY RADIANCE ON ZENITH ANGLE (IN SOLAR PLANE) 
FOR CLEAR SKY CONDITIONS. Wavelength = 0.45Mm; solar 
angle = 65.7° . 
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5 

Albedo Calculations 

In this chapter we shall calculate the albedo of the Earth-Atmosphere 
system using the radiative-transfer model and also develop mathematical 
relationships for energy conservation in realistic atmospheres. 

5.1 Albedo 

Albedo means the ratio of output power of a system to the input power 
and must therefore range in value from zero to one. In this report we 
consider three albedos and it is important that we define them carefully 
so that no confusion occurs. First, the single-scattering albedo, is 
the ratio of the volume scattering coefficient to the volume extinction 
coefficient and is a fundamental property of the medium at a given point. 

It can vary throughout the medium. Second, the surface albedo, p is the 
ratio of reflected power to incident power independent of an atmosphere. 
Third, we shall consider the total or planetary albedo, A of the Earth-Atmo- 
sphere system. The main object of the present investigation is to understand 
the effect which aerosols have on the albedo A. 

The satellite determination of the albedo or the radiant exitance of a 
planet is not an easy task. By definition, the albedo concept involves the 
total power integrated over a hemisphere and a satellite measurement can 
only be done at selected small areas over a planetary surface. Essentially, 
what is measured is the total spectral radiance. In this case, the radiance 
of an atmosphere increases with increasing turbidity if the surface albedo 
is low, whereas the radiance will decrease if the surface albedo is high. 

This effect is illustrated in Fig. 14 in which the various atmospheric 
states are designated by the total transmittance. 

The total flux incident upon a plane-parallel atmosphere is y o E Q 

where u is the cosine of the solar zenith angle and E is the extra- 
o o 

terrestrial solar flux at the top of the atmosphere on a plane perpendi- 
cular to the sun. The radiant exitance is E + (o) , the power density lost 
through the surface is (l-p)E (t q ) , and the power density lost by absorption 
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in the atmosphere is (l-oo )E„. In order that energy be conserved we must 

0 X/ 

have 


p E = E,(o) + (l-p)E (t ) + (l-aj ft )E ? 

O O "i ““ O OX/ 


(50) 


The albedo of the total Earth-Atmosphere system is then defined as 


A = E + (o)/v o E o 


(51) 


or 


A = 1 - 


(l-p)E_(T 0 ) 


y e 

o o 


- ( 1-00 ) 


o y E 

O 0 


(52) 


It is interesting to consider special cases of Eq. (52) . If the 
atmosphere absorbs no radiation, i.e. 00^ = 1 then, using the specific 
relations for fluxes given in Appendix I we have for a pure scattering 
atmosphere , 


V + (l-n)T 0 | 

i + 2(i-n) ( i-p)t q 

y o + d-n)x o ' 

ji + 2(i-n) ( i-p)t q 


(53) 


and, as the atmosphere becomes optically very thin, i.e., as + 0 we 
see from Eq. (53) that A p. 

Another interesting case is when the atmosphere becomes optically 
very thick, i.e. as T q > 00 . Here we have 


(i-n)w 

o 

l-nw o + V(l-oo 0 ) (l+oo 0 - 2n"w 0 ) 


(54) 
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If, in addition, 0) = 1, we see that A = 1 as it should be. Such an 
atmosphere behaves as a perfect reflector. On the other hand, if 0) = 0 
for this case we see that A = 0, a perfect absorber. A very simple re- 
lationship holds if we have a Rayleigh or an isotropic scattering law, 
i.e. f| = 1/2. In this case 


01 

A 2 (55) 

2-oi + 2\/l-oi ) 
o o 

Some of these relationships can be illustrated as in Fig. 15. We 
have calculated the total albedo as a function of atmospheric optical 
thickness and surface albedo for a pure scattering atmosphere. The reason 
for the minimum is that the atmosphere was gradually changed from a 
Rayleigh atmosphere with p = 1/2 to an atmosphere with a greater aerosol 
contribution in which p > 1/2. 

5.2 Energy Conservation Relations 

Equation (52) can be rewritten as follows: 

A + (l-p)F + (l-oi )F„ = 1 (56) 

— OX/ 

This conservation relation can be illustrated graphically by considering 
each component to be the length of a perpendicular in an equilateral 
triangle of unit altitude. Such a relationship is illustrated in Fig. 16. 
Any planetary state can be represented by a. point located somewhere 
within the triangle. A planetary system with low surface albedos and 
very little absorption will have points in the lower right corner of the 
diagram. A system with low surface albedo and strong absorption will 
have points along the left side of the diagram and systems with high 
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surface albedos and very little absorption will have points in the 

upper corner of the diagram. Hence, if a series of measurements are taken 

the energy loss of an atmosphere can be estimated. 

If we consider non-absorbing atmospheres (to = 1) and separate the 

total downward flux F into a direct solar component F, and a diffuse 

— d 

downward flux F then we again have a conservation relation, 

A + (l-p)F + (l-p)F = 1 (57) 

d — 

As in the other diagram this relation can also be represented by 
points within an equilateral triangle. It is illustrated in Fig.. 17 
wherein any non-absorbing planetary atmosphere system can be represented. 

Thin atmospheres with low surface albedos are represented by points in the 
lower right corner of the triangle; thick atmospheres with low surface 
albedos have points along the left side of the triangle, and systems with 
high surface albedos are represented by points in the upper corner of the 
triangle. Using the radiative-transf er formulas I have calculated these 
components for various planetary systems with optical thickness ranging 
from a Rayleigh optical thickness to multiples thereof as high as 100. 

These cases are illustrated in Fig. 18. 

5.3 Albedo-Radiance Relationship 

Using the atmospheric-radiative- transfer model we can calculate the 
total albedo as a function of the total radiance. Such a relationship 
is more nearly consistent with an actual measurement of radiation by a 
satellite. Figure 19 illustrates the dependence of total albedo on the 
total spectral radiance in the vertical direction at the top of the 
atmosphere. Thirty-six planetary systems are represented, i.e. six surface 
albedos and six atmospheric states beginning with a pure Rayleigh atmosphere. 
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FIGURE 17 

ENERGY CONSERVATION RELATION FOR 
A NON-ABSORBING PLANETARY 
ATMOSPHERE 

A <v Total Albedo 
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FIGURE 18 

RELATIONSHIP OF TOTAL ALBEDO, 
DOWNWARD DIRECT FLUX LOSS RATIO, 
AND DOWNWARD DIFFUSE FLUX LOSS 
RATIO FOR ATMOSPHERES WITH NO 
ABSORPTION. SOLAR ZENITH 
ANGLE = 30°; WAVELENGTH = 0.45 ym. 

i - RAYLEIGH ATMOSPHERE 
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The arrow on each line indicates the progression from a Rayleigh to an 
extremely dense atmosphere. It should be noted that the relationship 
is by no means a simple one; the complexities involve the competing effects 
of surface albedo and atmospheric turbidity. It does demonstrate, however, 
that given the surface albedo and the radiance the total albedo can be 
determined. 

Another interesting relationship is the connection between total albedo 
and sky radiance. This is especially useful if a satellite or aircraft 
measurement is not feasible. The dependence of total albedo on vertical 
sky radiance is illustrated in Fig. 20. 

5.4 Albedo-Irradiance Relationship 

One of the easiest measurements to perform is that of the total 
irradiance at the surface. A flat-plate device will receive radiation from 
the entire hemisphere and include the direct solar as well as the downward 
diffuse irradiance. Again, model calculations were made to show the re- 
lationship between total irradiance and total albedo. This is illustrated 
in Fig. 21. The Rayleigh atmosphere is designated .by the letter R. Thus, 
we see that for a given surface albedo and a measure of the total irradiance 
the total albedo can be determined. 

A similar relationship exists between the diffuse sky irradiance 
and total albedo. This effect is illustrated in Fig. 22. In this case 
and in the preceding case it should be noted that there is a linear 
relationship between the total albedo and irradiance. The diffuse component 
of the flux is more difficult to measure since it requires that the sun 
be blocked out of the hemispherical field of view. 

5.5 Total Albedo-Surface Albedo Relationship 

Finally, we should understand the relationship between the total albedo 
and the surface albedo. The connection between the two is depicted in 
Fig. 23 in which the straight line represents zero optical thickness or no 
atmosphere. It should be noted that as the optical thickness of the atmo- 
sphere increases the total albedo decreases for a given surface albedo. 
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FIGURE 22 ' 

TOTAL ALBEDO VS. DIFFUSE 
IRRADIANCE FOR OPTICAL 
THICKNESSES OF 1, 2, 5, ll 
50, 100 x RAYLEIGH VALUE. 

SOLAR ZENITH ANGLE - 30' 
WAVELENGTH - 0.45 um 


DIFFUSE IRRADIANCE (mW/cm -wa) 


NO AEROSOL ABSORPTION 
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This is due to the increasing degree of anisotropy of the scattering as 
more particulates are added to the atmosphere. Beyond an optical thickness 
of about 10 t however, the degree of anisotropy of the scattering seems to 
have less effect and the total albedo is then determined primarily by 
the ever-increasing atmosphere. 
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6 

Experimental Data 

In this chapter we shall analyze the actual data obtained by SKYLAB 
and by other auxiliary means and make comparisons with the theoretical 
models. 

6.1 SKYLAB Data 

In this section we shall examine the SKYLAB data. The optimum site 
conditions for this investigation would have been to have a well-documented 
data collection area in or around the target area. The target and back- 
ground should ideally be black so as to reduce the direct influence of 
the surface. Thus, the remote sensing equation would read 

L = L p (58) 

i.e. we have only the atmospheric path radiance to consider. Unfortunately, 
the real world does not have zero reflectance surfaces. During the first 
SKYLAB mission it was planned that data would be collected by the S-192 
multispectral scanner over Pisgah Crater in California, a region with a 
very low surface reflectance. Difficulties with the original SKYLAB 
precluded this site however and the White Sands, New Mexico, site was 
chosen instead. On another SKYLAB mission we were able to get multi- 
spectral data for the Lake Michigan site. 

6.1.1 White Sands, New Mexico 

I chose two specific sites in the White Sands area. One was in the 
southern part of the bright White Sands area near Lake Lucero. The other 
was in the southern-most point of the lava flow and is representative of 
a region with a low surface albedo. These will be referred to as White 
Sands Bright and White Sands Dark respectively. The sites were chosen by; 
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looking at SKYLAB photographic data in order to find a cloud-free area. 

From an account of the weather at the time SKYLAB flew over this particular 
site it was learned that some clouds were present, a fact which was not ob- 
vious upon detailed examination of the photographic data. It is easy to 
confuse bright sandy areas with clouds. 

6.1.2 Lake Michigan 

Another site which I chose was Lake Michigan under clear skies and 
in the deep-water area. The particular region chosen was about 10 miles 
west of St. Joseph, Michigan near the southern tip of the lake. At the 
time of the SKYLAB data take there were some clouds along the shoreline 
but the target area was clear. Besides the SKYLAB multispectral data we 
also had an aircraft overflight using the multispectral scanner aboard 
the ERIM-owned aircraft. In addition, the plane had an integrating 
nephelometer on board which sampled the air to collect data on the 
volume extinction coefficient. 

6.2 Model Comparisons 

We can now look at the absolute radiance values and compare the 
experimental results with our theoretical calculations. 

6.2.1 White Sands, New Mexico 

The specific White Sands National Monument consists of several 

hundred square miles of nearly pure gypsum (CaSO^ • ^H^O) and because of 

its extensive area and high reflectance it is clearly visible at spacecraft 

altitudes. Careful measurements have been made of the spectral reflectance 

L 36 1 

of specific regions in the White Sands area by Lindberg and Smith who 
are permanently employed in the site area. The spectra for a sample of 
white sand, playa crust, and basalt are illustrated in Fig. 24. These 
spectra were used in the model comparison. Table II lists the S-192 bands, 
the corresponding wavelengths of the band centers, the spectral reflectance 
of white sand, the SKYLAB radiances, the experimental standard deviation, 
and the theoretical value of radiance as calculated from the radiative- 
transfer model. The model atmosphere chosen was a pure Rayleigh atmosphere 
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TABLE II. COMPARISON OF SKYLAB DATA WITH MODEL EOR THE WHITE SANDS 
(Bright) AREA. 

Date: 14 June 1973 ; GMT: 14:44:46 ; Solar Zenith Angle: 32°28’54.5" ; Nadir View Angle: 5° 

Azimuth Angle: 47°6 , 36" ; Base Altitude: 1.2393 km 


Band 

Number 

Wavelength 

(Um) 

Surface 

Reflectance (%) 

SKYLAB Radiance 
(mW / cm 2 -sr-ym) 

Standard 

Deviation 

Model Radiance 
(mW / cm 2 -sr-ym) 

1 

0.4335 

55.0 

5.873 

0.3426 

23.84 

2 

0.4770 

60.0 

10.515 

0.4458 

32.23 

3 

0.5250 

62.5 

11.959 

0.8167 

30.23 

4 

0.5700 

67.5 

17.473 

1.648 

30.17 

5 

0.6265 

70.5 

13.894 

0.7837 

29.09 

6 

0.6940 

71.0 

10.968 

0.4972 

25.71 

7 

0.8300 

73.0 

9.341 

0.4957 

19.54 

8 

0.9900 

70.0 

5.087 

0.2262 

13.93 

9 

1.1100 

74.0 

2.942 

0.1169 

11.17 

10 

1.2150 

60.0 

2.957 

0.1380 

7.34 

11 

1.6400 

52.0 

1.160 

0.05514 

3.08 
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FIGURE 24. Diffuse reflectance spectra of reagent grade CaSO^ 2^0 along 
with typical examples of white gypsum sand, Lake Lucero playa crust, and 
basalt. The curves have been smoothed to eliminate spectrophotometer noise 
and baseline effects. 


which should give us the lowest radiance values. Instead, it should be 
noted that they are all consistently higher than the SKYLAB values by 
a factor of approximately three. This is illustrated in Fig. 25- How such 
a large discrepancy can exist is very puzzling. It is interesting that the 
ratio of theoretical-to-experimental values is nearly equal to pi. We 
took the theoretical radiances and divided them by pi and the experimental 
values. The spectrum is plotted in Fig. 26. Although there is a large 
fluctuation in the ratios the mean seems to be about 0.95. Could it be 
that the original calibration was off by a factor of pi? The experimental 
values were checked several times by individuals at ERIM who were thoroughly 
familiar with the SKYLAB multispectral calibration procedures and no dis- 
crepancy was indicated. I discussed the matter with Lawrence Korb now 
at the Goddard Spaceflight Center in Greenbelt, Md. He indicated to me 
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that there were probably calibration problems somewhere. Is there any 
way of resolving the matter? If we include absorption in the model 
calculations it reduces the theoretical values by only a very small amount, 
certainly not enough to bring theory into agreement with experiment. Only 
an excessively large optical thickness could bring the two values into 
agreement as is illustrated in Fig. 27 . It should be noted that by com- 
paring the SKYLAB radiances in Table II with model extrapolation in Fig. 27 
we would get optical thicknesses greater than 10, a completely unreasonable 
value. Is it possible that the model is incorrect? We have used this 
atmospheric-radiative-transfer model in conjunction with other experiments 
involving aircraft and LANDSAT data analysis and we have found no serious 

discrepancy. I have also compared the results of the model calculations with 

r 35 1 

those of Coulson et al. and the agreement is excellent. In addition, 

[37i 

Coulson and Jacobowitz calculated the radiance at the top of the atmo- 

sphere specifically for the White Sands area and got the same results that 
I did. Thus, there are three independent calculations of the spectral 
radiances which agree but which disagree with the SKYLAB values. 

I have discussed the lack of agreement with J. Lindberg and E. Williamson 
both of the Atmospheric Sciences Laboratory at White Sands, N.M. Lindberg 
said his measurements have indicated an imaginary part of the particulate re- 
fractive index to be about 0.01 in the visible and near-infrared part of the 
spectrum, a value which is still too small by itself to decrease the 
theoretical values significantly. It was learned however, that it did rain 
in the site area prior to the SKYLAB data take and this could have reduced 
the surface albedo. It is possible that the reflectance could therefore be 
lower than the dry values by perhaps a factor of two. Thus, we conclude 
that only if the reflectance were reduced by an unreasonably large factor 
and the atmosphere at the time of the SKYLAB data take were strongly 
absorbing and optically very thick could there be agreement. 
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FIGURE 27 

EXTRAPOLATION OF THEORETICAL 
CALCULATIONS TO LARGE OPTICAL 
THICKNESSES FOR THE 
WHITE SANDS BRIGHT AREA 
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We now turn to the analysis of SKYLAB data for the White Sands Dark 
area, i.e. the basalt region in the lava flow just north of the White Sands. 
Table III lists the comparisons between SKYLAB data and model calculations. 

A graphical display is depicted in Fig. 28. This time the model gives 
values about 1.5 times greater than the experimental radiance values and 
again it is difficult to see how they can be brought into agreement. It 
is true that wet basalt would have a reflectance smaller than the one used 
in the model calculations but not small enough to agree, with experiment. 
Also, the lack of agreement seems to be about the same across the entire 
spectral region. We conclude that the same conditions prevail for this 
site as for the previous one. 

6.2.2 Lake Michigan 

Estimates of the reflectance of water are more difficult to make 
than for dry land as a result of the large temporal variations which can 
exist. The data which we chose were taken from measurements of Santa 

r 88 1 

Monica by Wezernak L ' . A comparison of model calculations with SKYLAB 
data are listed in Table IV and a graph is illustrated in Fig. 29. Sur- 
prisingly, the agreement is much better than for the White Sands area 
although it should be noted that it is possible to make almost any 
theoretical calculation go through the experimental points. The experi- 
mental error is very large as can be seen in Fig. 29. 

The ERIM multispectral scanner was flown over the site at almost 
the time of SKYLAB overpass and data were collected at three altitudes; 
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TABLE III . COMPARISON OF SKYLAB DATA WITH MODEL FOR THE WHITE SANDS 
(Dark) AREA. 


Date: 14 June 1973 ; GMT: 14:44:42 ; Solar Zenith Angle: 32 o 20’7.25" ; Nadir View Angle: 5°32' ; 

Azimuth Angle: 24°52 , 9.623" ; Base Altitude: 1.280 km 

Band Wavelength Surface SKYLAB Radiance Standard Model Radiance 

Number (ym) Reflectance (%) ( mW / cm 2 -s r-ym) Deviation (mW/cm 2 -sr ym) 


1 

0.4335 

4.0 

2.847 

0.1959 

4.954 

2 

0.4770 

4.5 

3.574 

0.2190 

5.171 

3 

0.5250 

4.9 

2.939 

0.2885 

4.065 

4 

0.5700 

5.0 

2.071 

0.6202 

3.368 

5 

0.6265 

5.1 

1.656 

0.7079 

2.823 

6 

0.6940 

5.5 

1.091 

0.3422 

2.408 

7 

0.8300 

6.0 

1.375 

0.3518 

1.770 

8 

0.9900 

7.0 

0.927 

0.2636 

1.453 

9 

1.1100 

7.5 

0.562 

0.1372 

1.162 

10 

1.2150 

7.9 

0.572 

0.2015 

0.985 

11 

1.6400 

8.8 

0.303 

0.0811 

0.524 
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TABLE IV. COMPARISON OF SKYLAB DATA WITH MODEL FOR THE LAKE MICHIGAN AREA 


Date: 18 Sept 1973 ; GMT: 16:00:11 ; Solar Zenith Angle: 43°41 T 13.25 n ; Nadir View Angle: 5°32' ; 

n-an” . n 1-7-7 1,™ 


Azimuth Angle: 102°48’17. 

039" ; Base Altitude: 

0.177 km 



Band 

Number 

Wavelength 

(ym) 

Surface 

Reflectance (%) 

SKYLAB Radiance 
(mW / cm 2 -s r -ym) 

Standard 

Deviation 

Model Radiance 
(mW / cm 2 -s r-ym) 

1 

0.4335 

3.4 

3.107 

0.1888 

4.419 

2 

0.4770 

2.6 

4.630 

0.1894 

4.095 

3 

0.5250 

2.35 

3.849 

0.2326 

3.035 

4 

0.5700 

1.9 

5.488 

0.3448 

2.417 

5 

0.6265 

1.6 

1.252 

0.4033 

1.653 

6 

0.6940 

1.3 

0.5681 

0.2258 

0.971 

7 

0.8300 

1.0 

0.2979 

0.2557 

0.402 

8 

0.9900 

0.5 

0.1276 

0.2337 

0.149 

9 

1.1100 

0.4 

0.1095 

0.1309 

0.083 

10 

1.2150 

0.3 

0.0644 

0.2331 

0.049 

11 

1.6400 

0.2 

0.0197 

0.1088 

0.013 
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1000, 5000, 10,000 ft. These data with SKYLAB computations are plotted 
in Fig. 30, It is interesting that the spectra have almost the same 
spectral shape, They also become progressively larger in absolute radiance 
value as we go to higher altitudes, as should be. Unfortunately, one cannot 
say anything about the optical properties of the intervening atmosphere even 
though measurements were made of the volume extinction coefficient at the 
three aircraft altitudes. The reason is that the integrating nephelometer 
indicated a higher reading at 10,000 ft than it did at 5000 ft, a fact which 
is difficult to explain, Also, since there are no measurements at the sur- 
face or at altitudes greater than 10,000 ft there is no way that an optical 
thickness can be estimated. In the model calculations, we therefore assumed 
a pure Rayleigh atmosphere. Nevertheless, we also calculated the radiance 
assuming a small amount of absorption and the difference was very small. 
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Conclusions and Recommendations 

In this investigation we have developed new and simple procedures 
for calculating and interpreting the spectral albedo of our planet. The 
triangular plots for the energy conservation relations are especially in- 
teresting and this method of representing atmospheric states shows promise 
for future studies of the atmospheres of planetary bodies. 

Methods were also indicated whereby one can now perform simple 
measurements of various radiances and irradiances either at the top or 
bottom of the atmosphere and obtain knowledge of the total albedo. 

It is unfortunate that the SKYLAB experimental data were unreli- 
able. About the only thing which can be said is that spectrally these data 
did seem to match the model calculations reasonably well. As indicated 
in Chapter 6, the SKYLAB radiances were all about one-third of the model 
calculations for the White Sands Bright area and about two-thirds the model 
calculations for the White Sands Dark area. This is extremely difficult to 
explain other than to assume that SKYLAB calibration problems existed. In 
order to bring any model calculation into agreement with experiment it is 
necessary to make unreasonable assumptions regarding both the state of the 
atmosphere and the character of the Earth’s surface. A different problem 
exists in the interpretation of the SKYLAB data for the Lake Michigan site. 
Here the agreement seems to be good but the SKYLAB error is so large , 
especially for the longer wavelengths that one can fit almost any atmo- 
spheric state to the experimental data points. Nevertheless, the aircraft . 
and SKYLAB data do agree very well spectrally. 

In spite of the difficulties associated with the SKYLAB data, I do 
believe that this method of determining the spectral albedo of our planet 
is a reasonable one. Eor future investigations, I strongly recommend that 
the spacecraft sensors be calibrated in the laboratory and periodically in 
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orbit by comparing radiance values with measurements made at selected test 
sites. Careful weather records should be taken at the test sites and 
radiance and irradiance should be measured periodically to insure the 
environmental stability of the area. If ground measurements and space- 
craft measurements are taken at the same time for many selected test sites 
and these measurements are repeated many times then more reasonable esti- 
mates could be made of the spectral albedo of Earth, Such a determination 
of the albedo of the Earth-Atmosphere system is more than of academic 
interest. Indeed, it may determine the fate of civilization on our planet. 


I 
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Appendix I 

We present here the solutions for the diffuse upward flux E (x) , 
the diffuse downward flux E_(x) , and the total (direct solar plus diffuse) 
downward flux E (x) for a multiple scattering and absorbing medium. 

E + (X) 

E_(X) 

E (X) 


y e 

o o 

[ (1-nw ) sinh kx + ky cosh kx ] 

1 v 1 o o o o 

wy 

[a-ni 0 o ) - pw 0 (l-n)] sink 2kW 0 t 0 + W 0 cosh 2ky 0 T 0 


where 

r = 
a = 


r |(l-n)03 o sinh k(x Q -x) + A j\l-nuj o )sinh 2ky Q X + ky Q cosh 2ky Q xjj 
T j(l-ya) o )sinh k(x Q -x) + ky Q cosh k(X Q -x) + A(l-n)w o sinh 2 ky Q xJ 


-t/y. 


- y E e 
o o 


= T |(l-nw ) sinh k(x -x) + ky cosh(x -x) + A(l-n )(0 sinh 2ky x 
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Appendix II 

The solutions of the radiative-transfer equation for path radiance 
and sky radiance for a multiple-scattering and absorbing medium are as 
follows : 


L p (T,y s (j)) = 


jsii 


Ao) Q (l-n) jsinh k(T Q -T) + ky 


’ /y 

e - cosh 1c(t 


o- T) ] 


p(u,4>,p 0 »tt+4> 0 ) 


+ A j(l-na) o -k 2 P o y) sinh k(T Q -T)V p ( y , 9 ) 


+ Ak| (l-no) )y-y 

o o 


+ $|l + 2(l-co o )y 


-(x o -T)/y 

e - cosh k(T -t) 

o 


p(y,cf) J -y o ,(f) o ) 


l+o) -2rio) 

o o 


-(T “T)/U 

sinh 2ky T~ e sinh 2ky T 

o oo 


+ ky $ 

o 


1 + 2(1+0) ~2yo) )y 

o o 


/y 

cosh 2ky T-e cosh 2ky T 

o o o 


and 


L s (x,-y,cj)) = 


Ao) (l-p) -j sinh k(x -t) + ky cosh k(T -t) lp(-y,cj),-y ,ir+d> ) 
of o o l o o 

~T/\ir 

-Ao) 0 (1-T])e sinh kX Q + ky cosh kT I p (— y , 0 » 7T+ ^ ) 0 ) 

r ~ T/y 1 

+A (1-nw + k 2 y y) Sinh k(r -t) - e sinh kT p (~y ,<j>,-y^ ,<P) 

o o L o o j oo 

+Ak j^(l-r|o) o )y + -oil cosh k ( t q — t ) — £ cosh kxj p(-y,cj),-y o >(j) o ) 

+ $[l - 2(l-o) o )y][l + 0) Q - 2yo)J sinh 2-ky Q X 


+ky o $ jl - 2(l+o) o -2yo) o )yjJ^ cosh 2ky Q X-e ] 
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A 


to E 

o o 


4tt £ l- (ky) 2 J ^(1-T1W Q ) sinh kx^ + kU Q cosh kxj 


$ = 


2to Pky E 

o o o 


4tt jjL— (2TclJ o U) 2 J £(l-r|to o ) sinh kx^ + kh Q cosh kxj 


'ky cosh 2ky x + ri-pto - (l-p)rito 1 sinh 2ky x 
I o o o ^ o o 00 
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